Astronomy <fe Astrophysics manuscript no. aal7459-ll © ESO 2011 

October 3, 2011 



Letter to the Editor 

The Baade-Wesselink p-factor applicable to LMC Cepheids* 

N. Nardetto^, A. Fokin^'^, P. Fouque^, J. Storm^, W. Gieren^, G. Pietrzynski^'^, D. Mourard^, P. Kervella'^ 

^ Laboratoire Fizeau, UNS/OCA/CNRS UMR6525, Pare Valrose, 06108 Nice Cedex 2, France 

e-mail: Nicolas .NardettoOoca. eu 
^ Institute of Astronomy of the Russian Academy of Sciences, 48 Pjatnitskaya Str., Moscow 109017 Russia 
^ Observatoire Midi-Pyrenees, Laboratoire d'Astrophysique, UMR 5572, Universite Paul Sabatier - Toulouse 3, 14 

avenue Edouard Belin, 31400 Toulouse, France 
^ Leibniz-Institut fiir Astrophysik Potsdam (AIP), An der Sternwarte 16, D-14482 Potsdam, Germany 
^ Departamento de Astronomia, Universidad de Concepcion, Casilla 160-C, Concepcion, Chile 
j«— «, ' ® Warsaw University Observatory, Al. Ujazdowskie 4, 00-478, Warsaw, Poland 

^ I LESIA, Observatoire de Paris, CNRS UMR 8109, UPMC, Universite Paris Diderot, 5 place Jules Janssen, 92195 

, Meudon, France 

Ph. 

O ■ Received ... ; accepted ... 

o _ 

. Context. Recent observations of LMC Cepheids bring new constraints on the slope of the period-projection factor relation 

(hereafter Pp relation) that is currently used in the Baade-Wesselink (hereafter BW) method of distance determination. 
The discrepancy between observations and theoretical analysis is particularly significant for short period Cepheids 
Aims. We investigate three physical effects that might possibly explain this discrepancy: (1) the spectroscopic S/N that 
, is systematically lower for LMC Cepheids (around 10) compared to Galactic ones (up to 300), (2) the impact of the 

metallicity on the dynamical structure of LMC Cepheids, and (3) the combination of infrared photometry /interferometry 
with optical spectroscopy. 

Methods. To study the S/N we use a very simple toy model of Cepheids. The impact of metallicity on the projection 
^ , factor is based on the hydrodynamical model of 5 Cep already described in previous studies. This model is also used 

' to derive the position of the optical versus infrared photospheric layers. 

' Results. We find no significant effect of S/N, metallicity, and optical- versus-infrared observations on the Pp relation. 

I Conclusions. The Pp relation of Cepheids in the LMC does not differ from the Galactic relation. This allows its universal 

application to determine distances to extragalactic Cepheids via BW analysis. 

Key words. Stars: oscillations (including pulsations) 
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1. Introduction lation). Their relation is at 2 sigmas of the semi-theoretical 



ABSTRACT 



o 



, ^ , relation of Nardetto et al. (2007, 2009). For instance for 
, For decades the Cepheid stars have been used to calibrate ^ ^^p^ ^^^^ p = 1.41 ± 0.05, while the theoretical value 
the distance scale and the Hubble constant through their -^^^ ^ ± 0.05 (Nardetto et al. 2009). An effect of the 
O well-known period-lummosity {PL) relation (Riess et al. physical nature of LMC Cepheids on the slope of the Pp 
^ . 2009a, 2009b, 2011, and Freedman et al. 2010 for a review), relation is not excluded and should be investigated. This 
T-H . Recently, using the Baade-Wesselink (BW) method to de- question of the universality of the Pp relation is critical 
^ ; termme distances of Cepheids, Storm et al. (2011b, second because it is a common assumption in the BW analysis of 
• ^ , paper of the series) has found that the if-band PL relation extragalactic Cepheids 
■ is nearly universal and can be applied to any host galaxy 
H [ regardless of its metallicity. 
. - - 1 The projection factor is a key quantity of the BW meth- 
ods: it is used to convert the radial velocity into the pul- 
sation velocity of the star. In their first paper Storm et The spectroscopic S/N is systematically lower for LMC 
al. (2011a) directly constrained the period-projection fac- Cepheids (around 10) than Galactic ones (up to 300). In 
tor (Pp) relation from observations. The zero point is based Sect. 1, we investigate the impact of S/N on the radial 
on the HST trigonometric parallaxes of Galactic Cepheids velocity, hence on the projection factor. In Sect. 2, we de- 
(Benedict et al. 2007, Groenewegen 2007, Fouque et al. termine whether the projection factor should be corrected 
2007, Merand et al. 2004), and the slope is derived from BW when observing LMC Cepheids of lower metallicity. Then, 
distances of LMC Cepheids (all Cepheids in the LMC used when combining an optical spectroscopic determination of 
by Storm et al. can be assumed to be at the same distance, the projection factor with an angular diameter derived from 
leading to an extra constraint on the projection factor re- m/rared observations, a mismatch of the position of the pho- 
tosphere is not excluded. We investigate this possibility in 



* Based on observations made with ESO telescopes at the Silla Sect. 3. using new CRIRES spectroscopic observations and 
Paranal Observatory under programme ID 280.A-5018(A) hydrodynamical modelling. 
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Fig. 1. (a) The modelled random spectral line profile de- 
rived from our toy model (with T4ad = 20 km/s) together 
with the degraded profile (S/N of 10). The vertical lines 
correspond to the centroid (solid) and Gaussian (dashed) 
radial velocities. One hundred random spectral lines are 
calculated, (b) The relative uncertainties of the radial ve- 
locities as a function of the S/N (same legend), (c) The 
relative uncertainty on the bi-Gaussian line asymmetry as 
a function of the S/N. 



2. S/N and p-factor 

The amplitude of the radial velocity curve depends on the 
method used to derive the position of the spectral line pro- 
file. Two main methods are currently used: the Gaussian 
fit and the first moment (or centroid). Only the first mo- 
ment method is insensitive to the width of the spectral line 
(turbulence and/or stellar rotation) and is thus appropriate 
for comparing the dynamical structure of Cepheid's atmo- 
sphere (Burki et al. 1982, Nardetto et al. 2006). However, 
these two methods are not suitable for LMC Cepheids with 
spectra of very low S/N (around 10). The method cur- 
rently used for distant (and thus faint) Cepheids is to cross- 
correlate thousands of spectral lines, derive a residual spec- 
tral line and apply the Gaussian fit method to determine 
the radial velocity (Baranne et al. 1986). Compared to a 
single line measurement of the radial velocity (using the 
first moment method), the projection factor is actually re- 
duced (regardless of the period of the Cepheid) by about 
5% when using the cross-correlation method (Nardetto et 
al. 2009). The S/N might affect these different methods of 
radial velocity determination. 

We want to verify this assumption using a toy model 
already described in Nardetto et al. (2006). The modelled 
spectral line profile (Fig. la) is randomly degraded from an 
S/N of 300 to 10, by steps of 10. This is done one hundred 
times. For each S/N, we derive the mean radial velocity 
corresponding to the Gaussian and centroid methods (av- 
eraged over the 100 profiles), together with their respective 
mean uncertainties. A first result is that the mean radial 
velocity derived from both methods (Gaussian fit and first 



moment) is insensitive to the S/N. This is a good indication 
that the cross-correlation method is probably not biased by 
the S/N. 

Additional results have to be mentioned. The first mo- 
ment method can lead to statistical uncertainties as large 
as 20% when the S/N is lower around 10 (Fig. lb). The 
Gaussian fit method is more robust with a few percent 
precision for an S/N of 10. One should, however, consider 
that these results depend mainly on the radial velocity con- 
sidered (i.e. the pulsation phase considered) and the spec- 
tral line depth. When using the cross-correlation method, 
the relative uncertainty on the radial velocity can drop to 
less than 1%. We also take the opportunity to show that 
line asymmetry (derived using the bi-Gaussian method de- 
scribed in Nardetto et al. 2006) is highly sensitive to the 
S/N (Fig. Ic). While the mean value of the spectral asym- 
metry is stable with the S /N (as found for the radial veloc- 
ities), an S/N of 100 is required to derive a line asymme- 
try at a 10% relative precision level. Studying the dynam- 
ical structure of pulsating stars atmospheres thus requires 
a minimum S/N of 100. The conclusion of this section is 
that the S/N should not affect the period-projection factor 
relation of Galactic versus LMC Cepheids. 

3. Metallicity and p-factor 
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Fig. 2. The decomposition of the projection factor. The 
filled and open diamonds correspond to the Milky Way and 
LMC models of S Cep. For a given model, the lower (resp. 
higher) values are the amplitudes of the radial (resp. pulsa- 
tion) velocity curves for each spectral line considered. The 
red and blue arrows indicate the decomposition of the pro- 
jection factor for the MW and LMC respectively: t (Po), 
^ (/grad), I (/og) and \ {p = Po/grad/og)- For detailed 
explanations, refer to Nardetto et al. (2007). 

The metallicity of Cepheids has an impact on their k- 
mechanism that drives the pulsation through an opacity 
process. As a consequence, the amplitude of the pulsation 
velocity of the star is affected by the metallicity: The lower 
the metallicity, the larger the amplitude of the pulsation ve- 
locity (Klagyivik & Szabados 2007; Klagyivik & Szabados 
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2009). But in principle, if the amplitude of the pulsation 
velocity is higher, the radial velocity should also be larger 
in the same proportion, and the projection factor p = 
should not change. The key-point is then to know whether 
the velocity gradient in the atmosphere of the star changes 
with the metallicity in such a way that the projection factor 
for a given spectral line will change. To answer this ques- 
tion, we consider the hydrodynamical model of 6 Cep in 
the Milky Way (Nardetto et al. 2004), based on the code 
by Fokin et al. (2004), and compare it to a consistent LMC 
model of the same star. The fundamental parameters of the 
LMC model of 5 Cop arc determined as follows. Following 
Luck et al. (1998) and Romaniello et al. (2008), we first as- 
sume Z = 0.008 (very close to their respective values of 0.01 
and 0.009). Wc then consider a helium-to-metal enrichment 
ratio ^ of 1.5 (consistent with the value of 2.1±0.9 derived 
by Casagrande et al. 2007), which leads to X ~ 0.730 and 
Y = 0.262. Then, we assume the same mass for S Cep in the 
Milky Way and the LMC (i.e. M = 4.8M©) and we simply 
use the mass-luminosity relation of Chiosi et al. (1993, Eq. 
25) to derive the luminosity of S Cep in the LMC. Wc ob- 
tain L = 229OL0. Finally, using the Milky Way and LMC 
analytical relations of the red edge instability strip by Bono 
et al. (2000, Tab. 5), we find an expected effective temper- 
ature of about Tcff = GlOO-ft'. All parameters used in the 
code arc presented in Table 1 for the Galactic and LMC 
models of d Cep. We find a pulsation period of P = 5.38d 
(in LMC), which is consistent with theoretical results by 
Bono et al. (2000, Table 6). 

For both models (with Galactic and LMC metallicities), 
we calculated the spectral line profiles for the 17 metallic 
lines presented in Nardetto et al. (2007, hereafter Paper I) 
using radiative transfer equations in a moving atmosphere. 
The calculation is made at 40 equidistant pulsation phases. 
We finally plotted the amplitude of the radial (using the 
first moment method) and pulsation velocity curves as a 
function of the line depth at minimum radius (see Fig. 2) in 
order to measure the velocity gradient and perform the de- 
composition of the projection factor as explained in Paper I. 
There are three subconcepts involved in the decomposition 
of the projection factor: 1/ the geometric projection fac- 
tor mainly related to the limb-darkening of the star (po), 
2/ the correction due to the velocity gradient within the at- 
mosphere (/grad) and 3/ the correction due to the relative 
motion between the optical and gas layers corresponding 
to the photosphere (/og); for a detailed definition of these 
two layers, refer to Eqs. 1 and 3 of Nardetto et al. (2004). 
This decomposition of the projection factor for the Milky 
Way and the LMC models of 5 Cep are indicated in Fig. 2, 
and the corresponding values can be found in Table 1. To 
compare the MW and LMC projection factors properly, we 
have to consider two spectral lines of similar depth (around 
0.1): Fel 4896.439A for the MW and Fel 5373.709A for the 
LMC. These choices only have an impact on /grad- 

We find that the projection factor is about the same for 
the Galactic and LMC models (at a 1.5% level). This ba- 
sically means that even if the amplitude of the radial and 
pulsation velocity curves are larger for the LMC metallic- 
ity model, the structure of the atmosphere is relatively un- 
changed: the velocity gradient is larger (in absolute terms), 
but it is the same relative to the amplitude of the veloci- 
ties. As a conclusion, the p-factor value appropriate for a 
Cepheid of a given period should not depend on its metal- 



Table 1. Hydrodynamical models and corresponding p- 
factors 



parameters 


MW 


LMC 


L [Lq\ 


1995 


2290 


M [Mq] 


4.8 


4.8 


Teff [K] 


5877 


6100 


X 


0.700 


0.730 


Y 


0.280 


0.262 


Z 


0.020 


0.008 


P [days] 


5.41 


5.38 


Po 


1.390 


1.390 


/grad 


0.990 


0.975 


line for ./grad [A] 


Fcl 4896 


Fel 5373 


/og 


0.963 


0.963 


P = Po/grad/og 


1.325 


1.305 



licity, and the large discrepancy obtained for the Pp re- 
lations (observation versus theory), especially for short- 
period Cepheids (for 6 Cep, p = 1.41 [observations] versus 
p = 1.25[theory]), is not explained by a difference of metal- 
licity. 

4. Optical p-factor versus infrared diameters 

When applying the BW method, infrared photome- 
try/interferometry is often combined with optical spec- 
troscopy to derive the distance of the star. In principle, 
the radial velocity curve derived from optical spectroscopy 
has to be combined with an optical projection factor to 
obtain the pulsation velocity associated to the optical pho- 
tospheric layer. As a consequence, to get an accurate deter- 
mination of the distance, the limb-darkened angular diam- 
eter derived from the infrared surface brightness method or 
from infrared long-baseline interferometry must correspond 
to the same optical photospheric layer. In other words, an 
implicit assumption in the BW method is to consider that 
the photospheric layer in the optical and in the infrared are 
associated with the same material in the star. However, the 
larger the wavelength, the lower the position of the photo- 
sphere in the star (see theoretical results by Sasselov et al. 
1990, Fig. 21 and 22). This effect seems to be stronger for 
long-period Cepheids (their Fig. 21), than for short-period 
Cepheids (their Fig. 22). Verifying these theoretical results 
with observations is extremely difficult. The best is prob- 
ably to use interferometry. From the spectrocopic point of 
view, one can compare the optical and infrared radial ve- 
locity curves of two spectral lines of similar depths (using 
the first moment method). However, the velocity ratio ob- 
tained has to be corrected from the velocity gradient and 
from optical-versus-gas layer effect (/grad and fog, respec- 
tively). For instance, it is probably not possible to perform 
this analysis for atypical Cepheids, like X Sgr presenting a 
double shock wave in the atmosphere (Mathias et al. 2006), 
as done by Sasselov et al. (1990). 

We present infrared CRIRES observations of £ Car ob- 
tained in 2008-2009 at four different pulsation phases (S/N 
of about 300 in the continuum). We find three unblended 
spectral lines, Br7 21661. 2A and two metallic lines: Fell 
21266.1A and Cal 22090.0A. In Fig. 3, we present the Cal 
spectral line flux as a function of time. The first moment 
radial velocity and depth (as a percentage of continuum) 
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are indicated in Fig. 4, together with HARPS observations 
of a spectral hne of similar depth. The dashed line in Fig. 4 
(lower panel) is the optical HARPS radial velocity curve 
multiplied by 1.30, showing that the infrared radial velocity 
curve has an amplitude that is about 30% larger than in the 
optical. The conclusion is that, even if these two metallic 
lines in the optical and in the infrared have a similar depth 
(which means in principle that their line-forming regions 
have a similar optical depths, and that they form approxi- 
matively at the same position in the atmosphere compared 
to the photosphere), the infrared amplitude of the radial 
velocity curve is significantly larger than the optical one. 
This suggests (as explained above) that the photospheric 
layer is probably forming higher in the infrared than in the 
optical (which seems to be in contradiction with Sasselov 
et al. 1990). More data are required to confirm this result. 

Wc address this qiicstion theoretically using our hydro- 
dynamical models of 6 Cep and I Car. The model of I Car 
is described in Paper I. We consider the solar metallicity. 
Concerning 5 Cep, the optical and infrared photospheric 
layers (tc = 1) form at the same layer in the atmosphere, 
whereas for I Car the infrared photosphere forms 2% higher 
than the optical one. This implies a 2% increase only in the 
projection factor, which seems to be in contradiction with 
Sasselov et al. (1990) and also CRIRES data. This delicate 
question will be studied in more detail later, using a larger 
sample of infrared (rcsp. optical) spectroscopic (resp. inter- 
ferometric) data. In either case, this effect (optical versus 
infrared) should affect the Milky Way and LMC Cepheids 
in the same way. 
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Fig. 3. The CRIRES spectral line profiles (Cal 22090.0A) 
of i Car as a function of the pulsation phase. 

5. Conclusions 

Wc find that the S/N, metallicity, and optical- vcrsus- 
infrarcid observations (but this last point has to be con- 
firmed) cannot explain the discrepancy between the theo- 
retical and empirical period-projection factor relation based 
on LMC observations. Other possibilities can be consid- 
ered, such as limb darkening for instance. Indeed, following 
the relation by Storm et al. (2011a), the very short-period 
Cepheids should have a projection factor close to 1.45-1.5. 
This means a limb-darkening close to zero (uniform disk). 
By constraining this limb darkening using interfcrometry 
(for instance VEGA/CHARA in optical, Mourard et al. 
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Fig. 4. The line depth (up) and the first moment radial ve- 
locity (down) corresponding to the infrared CRIRES spec- 
tral line of t Car (crosses) is compared to the optical 4896A 
HARPS spectral line of similar depth (solid line). The 
meaning of the dashed line is explained in the text. 

2009), one can derive a geometric projection factor that 
should help in resolving the discrepancy. However, our re- 
sults seem to indicate that the Pp relation is universal. 
This is extremely precious when applying the BW method 
to extragalactic Cepheids. 
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